Disks on a Sphere and two-dimensional Glasses 

Werner Krauth * 
CNRS-Laboratoire de Physique Statistique 
Ecole Normale Superieure, 
24, rue Lhomond, 75231 Paris Cedex 05, France 

February 1, 2008 



Abstract 

Talk given at the conference on 'Inhomogeneous Random Systems' 
at the University of Cergy-Pontoise, France, (23 January 2001). 

I describe the classic circle-packing problem on a sphere, and the ana- 
lytic and numerical approaches that have been used to study it. I then 
present a very simple Markov-chain Monte Carlo algorithm, which suc- 
ceeds in finding the best solutions known today. The behavior of the 
algorithm is put into the context of the statistical physics of glasses. 
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What is the minimal radius -Rtv of a sphere ahowing N non-overlapping 
circles of radius r = ^ to be drawn on its surface? In a classic paper, Schiitte 
and van der Waerden [Q] reduced this problem to a study of planar graphs: 
As the sphere has minimal radius, many circles are blocked, i.e. they cannot 
move. These circles make up the M < N vertices of the graph. If two circles 
are in (blocking) contact, they are connected by an edge. 

Schiitte and van der Waerden considered all graphs with M < 12 ver- 
tices. Some of these graphs can correspond to a blocked configuration of 
circles on a sphere of radius R, with R = miugraphs R- They were able to 
determine Rn for N < 12. For this latter case N = 12, e.g., the centers of 
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Figure 1: For = 13 spheres, 1 < = as proven by Schiitte and van 
der Waerden. This means that 13 identical spheres cannot be placed on top 
of another one, of same size. 



the circles are located at the vertices of an icosahedron of side length 1, and 
lie on a sphere of radius 

= = 0.951... (1) 

Schiitte and van der Waerden Q were later able to prove that R > 1 for 
> 12. This result (c/ figure |^) implies that 13 spheres of radius 1 cannot 
be placed on top of another sphere of the same size, and solves a packing 
problem discussed at least since the time of Newton. The best currently 
known packing of spheres for = 13 concerns spheres of radius 1 on top of 
another sphere of radius 1.0911 .... The corresponding circle arrangement is 
shown in figure |2|, and is almost certainly optimal. Today, it should be easy 
to enumerate planar graphs on a computer, and to prove this conjecture. 




Figure 2: Best known arrangement of 13 circles of radius ^ (on the surface 
of a sphere of radius 1.04557...). As a consequence, 13 spheres of radius 1 
can be placed on the surface of a central sphere of radius r > 1.0911 .... 

The problem of packing disks on the surface of a sphere has appeared 
repeatedly in the literature, as there are important applications in mathe- 
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matics and the natural sciences. Hundreds of papers have been written on 
this subject (for a partial list of references cf e.g. |^]), most of them improv- 
ing upper bounds on i?jv in the range of 13 < ~ 100. Only for N = 24: 
has a new optimal solution been proven Q. 

The problem of N non-overlapping disks on a minimal sphere can be 
reformulated as a minimization problem of the potential energy of N parti- 
cles i = 1, . . . ,N on the surface of the sphere of radius 1 with a two-body 
interaction ^ 

N-l N , ^ p 

i = l j=i + l V 

where r^j is the Euclidean distance between circles i and j. Therefore, cus- 
tomary minimization programs (such as the Newton-Raphson method) can 
be applied to find local - and hopefully global - minima for ever increasing 
values of p. In our first units (circles of radius ^ on a sphere of radius R), 
the minimal radius R]\f is given by 

Rn = (3) 

miiiij Tij 

in the limit p — > oo. 

Current records were set by Kottwitz |^], who used values of p in the 
astonishing range of p = 80 to p = 1 310 720. Running such a sophisticated 
Newton-Raphson program is a delicate task as one has to avoid getting 
trapped in local minima of the potential energy. Nevertheless, the strategy 
proved to be successful, as a large number of new values for were found, 
and as these values have so far stood the test of time. 

Recently, I have reconsidered the non-overlapping circle problem from 
the view point of Markov-chain Monte Carlo algorithms. Instead of search- 
ing for the minimum of a regularized energy, I move circles during a large 
number of time steps t = 1, . . . ,Tmax- At time t, a single, randomly cho- 
sen circle i is moved a tiny bit into a random direction. Specifically, let 
us denote the position of i at time t hy fi, and let the displacement vector 
be 6r = {6x, 5y, 6z) with Gaussian random numbers 6x etc. that have zero 
mean and variance 1. Then the following move is considered 

- = ^±|-i? (4) 

lU + dr\ 

If the move can be performed without violating the constraints, it is ac- 
cepted, and rejected otherwise. This procedure implements a standard lo- 
cal Metropolis algorithm for a constant stationary probability distribution: 
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Any non-overlapping configuration of circles should be visited with the same 
probability. 

In addition to performing the Metropolis algorithm at fixed R, I also 
implemented the "simulated annealing" procedure Every so often, the 
radius of the sphere is reduced by a very small amount, if this is possible 
without introducing overlaps between particles . 

I have found that this simple algorithm reproduces the best results known 
up to date (c/ |Q) for N = 5 ... 90, typically within a few minutes of com- 
puter time [0 . 

This observation may be of practical interest, as the Monte Carlo al- 
gorithm is much simpler than the previous approaches. However, it also 
underlines the quality of the bounds obtained by Kottwitz, which I have 
been able to reach, but not to improve on. 

In physical terms, the graph giving Rj\f is the ground state of the system 
(which does not have to be unique: for N = 15, e.g., two best solutions are 
known |^]). Following Stillinger and Weber |jll|, all the other graphs that 
can be put on the surface of a sphere (with larger radii R > Rn) may be 
called 'inherent structures'. These are the configurations that the simulated 
annealing algorithm may get trapped in. If the Monte Carlo algorithm is 
trapped in an inherent structure at radius R, it could by continuity also get 
trapped by the same inherent structure at a slightly larger radius R' ~ R. 
This means that, in the strict sense, the Monte Carlo algorithm is non- 
ergodic at least for all R' satisfying 

max R> R' > Rn. (5) 

graphs 

In this and similar systems, very little is known rigorously about the number 
of inherent structures as a function of radius (their density of states), and 
enumeration methods at small N should be very helpful. For physical appli- 
cations, the thermodynamic limit ^ oo is of prime interest. In this limit, 
{R' — R)/R certainly goes to zero. This means that strict ergodicity (on 
timescales that diverge with N) is reinstalled in the thermodynamic limit. 
Precise understanding of this complicated thermodynamic limit is lacking. 

If, in contrast, the simulated annealing algorithm got trapped in a non- 
optimal state (with R clearly larger than Rn) with probability close to one, 
we would say that the 'physical' system of disks on the surface of a sphere 
had a glassy phase. Such a phase seems to be absent for monodisperse circles. 
Our findings agree very well with what has been found in monodisperse hard 
disk and hard sphere systems on two- and three-dimensional tori [ p!2| , 
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In contrast, there is strong numerical evidence that polydispeise hard 
disks can be glassy: in a system of hard disks of varying radii on a torus 
[p^ , a comparable local Monte Carlo algorithm always gets trapped at finite 
density. Again, it is fundamentally unclear how this observation carries over 
into the thermodynamic limit. Empirically, diffusion times of particles seem 
to grow without bounds at a well-defined, finite density. 

In these two-dimensional systems, it has been found fl^ that other 
Monte Carlo dynamics do not suffer the same slow-down as the local al- 
gorithm. These methods also provide strong indications that the glass tran- 
sition in these two-dimensional models is a purely dynamic phenomenon, 
which does not manifest itself in equilibrium averages over the stationary 
probability distribution. 

Acknowledgements: I would like to thank L. Santen and S. Tanase-Nicola 
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